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ABSTRACT 

This paper investigates multiwavelength retrievals of median equivolumetric drop diameter 
D o suitable for drizzle and light rain, through collocated 355-/527-nm Micropulse Lidar Network 
(MPLNET) observations collected during precipitation occurring 9 May 2012 at the Goddard Space 
Flight Center (GSFC) project site. By applying a previously developed retrieval technique for infrared 
bands, the method exploits the differential backscatter by liquid water at 355 and 527 nm for water drops 
larger than =50 gm. In the absence of molecular and aerosol scattering and neglecting any transmission 
losses, the ratio of the backscattering profiles at the two wavelengths (355 and 527 nm), measured from 
light rain below the cloud melting layer, can be described as a color ratio, which is directly related to D 0 . 
The uncertainty associated with this method is related to the unknown shape of the drop size spectrum 
and to the measurement error. Molecular and aerosol scattering contributions and relative transmission 
losses due to the various atmospheric constituents should be evaluated to derive D 0 from the observed 
color ratio profiles. This process is responsible for increasing the uncertainty in the retrieval. Multiple 
scattering, especially for UV lidar, is another source of error, but it exhibits lower overall uncertainty 
with respect to other identified error sources. It is found that the total error upper limit on D 0 ap- 
proaches 50%. The impact of this retrieval for long-term MPLNET monitoring and its global data 
archive is discussed. 


1. Introduction 

The primary objective of the National Aeronautics 
and Space Administration (NASA)’s Micropulse Lidar 
Network (MPLNET; Welton et al. 2001), a member of 
the Global Atmospheric Watch (GAW) Aerosol Lidar 
Observation Network (GALION; Hoff and Boesenberg 
2008), is the depiction of cloud and aerosol particle dis- 
tributions, utilizing a global and federated network of 
ground-based lidar instruments. Enhancing the project, 
efforts are presently ongoing toward the integration of 
ultraviolet (UV; 355 nm), eye-safe, and autonomous lidar 
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produced by Leosphere (Lolli et al. 2011) to complement 
the historical array of visible (VIS; 523, 527, and 532 nm) 
and eye-safe micropulse lidar instruments (Spinhirne 
1993; Spinhirne et al. 1995; Campbell et al. 2002). The 
integration of these instruments would facilitate the 
emergence of a denser global network. When possible, 
collocated UV and VIS lidar profiling would further 
provide multispectral information on aerosol particle 
and cloud scattering. 

Covering a large fraction of the planet surface at any 
given time, stratocumulus and stratus clouds play an im- 
portant role in planetary radiation budget (e.g., Bennartz 
et al. 2013). These clouds are important because they 
redistribute the column water liquid content due to the 
transformation of water falling from clouds into vapor 
when precipitating (Harrison et al. 1990; O’Connor et al. 
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2005). Numerical models require an accurate represen- 
tation of the microphysical properties of these boundary 
layer clouds in order to explicitly resolve their impact on 
climate (Slingo and Slingo 1991; O’Connor et al. 2005). 
The purpose of this paper is to describe and evaluate 
a method for measuring light rain and drizzle drop me- 
dian equivolumetric size diameter D 0 from collocated 
MPLNET UV and VIS lidar profiles by taking advan- 
tage of the differential backscatter during precipitation 
at the two wavelengths, respectively. Raindrop size is 
of particular importance for controlling the rate at 
which they deposit through the atmosphere and how 
quickly they evaporate below cloud base (Westbrook 
et al. 2010). 

Though atmospheric lidar profiling is attenuation lim- 
ited under many common aerosol- and cloud-observing 
conditions (Sassen and Cho 1992), a scenario exacerbated 
by liquid water droplet presence given the relatively high 
scattering cross sections and droplet concentrations 
(e.g., Sassen 1978), such data have been effectively used 
for investigating melting-layer thermodynamics and the 
scattering properties of precipitating low-level liquid 
water clouds (e.g., Sassen et al. 2005, 2007; Campbell 
and Shiobara 2008; Di Girolamo et al. 2012). Our goal is 
to transpose a two-color lidar-based droplet measurement 
technique designed for infrared bands (IR; Westbrook 
et al. 2010) to the UV-VIS region and to evaluate per- 
formance so as to extend, when possible, MPLNET 
(http://mplnet.gsfc.nasa.gov) capabilities for studying 
stratocumulus and stratus cloud precipitation size 
distributions. 


shape of the distribution (for p = 0, the distribution 
reduces to a simple inverse-exponential distribution). 
This distribution is preferred to log normal when 
precipitation originates from the melting layer (Willis 
1984). The total lidar backscatter for a distribution of 
water drops (range independent here), neglecting mul- 
tiple scattering [a reasonable assumption for MPL in- 
struments, given their relatively narrow field of view 
(Spinhirne 1993), though less so for Leosphere, as dis- 
cussed further below], is then 


fKN Q ,D Q ,r) 


1 

477 
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where the backscattering cross section cr bk is calculated 
using the Mie-based method described by Wolf and 
Voshchinnikov (2004) that is optimized for relatively 
large particles with respect to the profiling wavelength 
and with a complex index of refraction for liquid water 
droplets of [1.35 + /(2.4 X 10 -9 )] and [1.33 + i(1.6 X 
10 -9 )] at 355 and 527 nm, respectively (Lynch et and 
Livingston 2001). The values of o-bk(-D) are integrated 
over the size distribution [Eq. (2)] for D and D 0 
ranging in the interval 50-1000 pm to compute [3 at both 
wavelengths. 

Shown in Eig. 1 (top) is the backscattering efficiency 
(or backscattering gain; unitless), Q bkg , solved as 
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2. Methods 

The theoretical radar backscattering cross section can 
be computed from Mie theory (Bohren and Huffman 
1983) for homogeneous water spheres. This solution 
is valid when the raindrop diameter does not exceed 
1.0 mm (Beard 1976) — a value beyond which a raindrop 
should be otherwise considered oblate. To calculate the 
theoretical unattenuated lidar backscatter coefficient, 
precipitation is represented with a gamma-based size 
distribution (O’Connor et al. 2005; Westbrook et al. 
2010) for the number concentration of drops d N with 
diameters between D and D + d D as 


for both 355 and 527 nm. Note that <2 bkg , computed over 
a size distribution, can differ greatly at the two wave- 
lengths for a given median equivolumetric diameter D 0 
(Fig. 1, bottom). Shown in Fig. 2 are corresponding ex- 
tinction efficiency parameters Q cxl [Eq. (4)] (top) and 
(7 ext computed over a size distribution (at 355 and 527 nm, 
respectively; bottom) with o- ex t the extinction cross sec- 
tion calculated using the same Mie-based method for the 
backscattering cross section. This term is solved as 

2ext(£) = T^- (4) 

4^ 


-J±(N 0 ,D,D 0 ,») = N 0 (£-) exp 


D 
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where N 0 controls the total concentration of drops for 
a given D 0 and where p, is the median equivolumetric 
diameter and the dispersion parameter that controls the 


In contrast to backscatter, o- ex t is mostly insensitive to 
the raindrop diameter, reaching approximately 2, the 
geometric optical limit, for D > 50 pm. Then, through 
the differential backscattering efficiency at the two wave- 
lengths, it is possible to retrieve the median equivolumetric 
diameter of the precipitation, defining a backscatter 
color ratio (C R , dB) as 
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Fig. 1. Lidar Q bkg as defined in Eq. (3), (top) derived from Mie theory for spherical water 
drops in l-/rm droplet size intervals and (bottom) (2bkg computed over a size distribution de- 
fined in Eq. (1) with fj. = 2. 


C R (D 0 ,^) = mogJj^A, (5) 

\ 1 rain / 

where f3 355 and yS 527 are the backscattering co- 
efficients due to the rain calculated with Eq. (2) at 355 
and 527 nm. As the droplet concentration does not 
change with the wavelength, C R only depends on the 
droplet diameter and not on N 0 . Interpreting C R in 
terms of D 0 requires some information on the shape of 
distribution and therefore is sensitive to the value 
chosen for /x. In Fig. 3, we show profiles for C R solved at 
lx = 0, 2, 4, . . . , 10. 

Range-dependent normalized attenuated total back- 
scatter coefficient [cf. Eq. (12)] measured by the lidar 
during precipitation can be substituted into Eq. (5) to 
estimate C R (r) at any range gate r. With C R solved, D 0 (r) 
can be estimated. A theoretical C R function is ultimately 
two dimensional, and thus it can be solved in lookup 
database form (depicted visually in Fig. 3) and applied as 

C R =f(P Q ,p). (6) 

theor 

For any measured C R (r), D 0 is then retrieved, performing 
a bivariate analysis minimizing the absolute difference 
between measured and theoretical C R over D 0 and the 


dispersion parameter. This bivariate analysis avoids the 
assumption of the dispersion parameter /x, which is re- 
lated to the precipitation type: 

D 0 (r) = f~ l [C R ( r),/x \ . (7) 

meas 

Welton and Campbell (2002), Campbell et al. (2003), 
and Lolli et al. (2011) describe the calibration procedures 
necessary for retrieving attenuated and corrected back- 
scatter coefficients for MPLNET lidar, including the UV 
instrument, using an independent measurement of total 
column aerosol transmission from a nearby sun photom- 
eter. During precipitation episodes and cloud presence, 
however, standard techniques are not possible. Instead, 
the two signal profiles are normalized to corresponding 
attenuated molecular scattering profiles, as solved in a 
precipitation and aerosol particle scattering-free portion 
of the profile. Assuming that all instrument calibration 
terms have been accounted for in the processed signal 
P{r), the total backscatter coefficient j S Atot (r) relates as 

P A (r)r 2 = C A 0 A (r)/3 A ( r)[T \ (r)] 2 , (8) 

/l /l /v /l tot tot 

where C A is the linear system calibration coefficient and 
[r Atot (r)] 2 is the total two-way signal transmission. In the 
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Fig. 2. Corresponding with Fig. 1 Q ext as defined in Eq. (4) (top) derived from Mie theory for 
spherical water drops in l-/rm droplet size intervals and (bottom) Q cxt computed over a size 
distribution defined in Eq. (1) with /r = 2. 


equation, 0 A (r) is the overlap function, and it equals 
unity when the field of views of the telescope and laser 
fully overlap. All terms are with respect to system wave- 
length. A correction for overlap is needed for MPL in- 
struments (Campbell et al. 2002), as the range at which 
the full overlap range is reached is much higher than the 
UV lidar (~6.0 vs 0.35 km; Lolli et al. 2008). 


Similar to Eq. (7) of Campbell et al. (2008), however, 
by normalizing the signal to that of attenuated molecu- 
lar scattering within a “clear sky” section of the profile 
(i.e., where total backscatter is exclusively due to mo- 
lecular scattering), C A in Eq. (8) includes the aerosol 
transmissivity up to the reference zone. It is possible to 
define an intermediate value C* such that 



Fig. 3. The 355-/527-nm lidar color ratio, defined in Eq. (5), solved as a function of D 0 . 
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Ct = C x [T x (r ief )f[T x (r ref )] * 2 =- 


P A 0ref >rei 
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Each signal transmission term is distinguished here as 
relating either to molecular, rain, or aerosol particle 
components, and r rc( refers to a rain and aerosol particle 
scattering-free “clear sky” reference range. 


At any point in the profile, the total backscatter co- 
efficient can be estimated from Eqs. (8) and (9) as 


P^V[T k (r ref )] 2 [T A (r f )] 

ram aer 

C A *[r A (r)] 2 [T A (r)] 2 [r A (r)] 2 ’ 

rain aer mol 


and then the measured C R after the lidar signal cali- 
bration is 


C*(r) = 101og 10 


(Vh55 t SWiss. (r)] 2 [T 355 (r)] 2 [r 355 (r)] 2 }{[E 527 (r ref )] 2 [T 527 (r ref )] 2 }\ 
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CR retrievals and relative measurement uncertainty 


To retrieve D 0 from Eqs. (5) and (11), the transmission 
due to aerosol particles and rain and molecular and aerosol 
backscattering at two different wavelengths must be 
evaluated. In Eq. (11), for D 0 > 50 /nil. attenuation due to 
precipitation or [7) rain (r rc f)] 2 differs by less than 2% be- 
tween the two wavelengths. The transmission term due to 
the rain then can be simplified. Further, we assume that the 
aerosol particle mass loading is constant between the time 
of signal calibration and during the rainfall event. The 
aerosol transmission is evaluated by retrieving the back- 
scattering coefficient profile at calibration time, inverting 
the signal with a Fernald method (Femald 1984), and as- 
suming a fixed value for the ratio of extinction and back- 
scatter cross sections [i.e., the lidar ratio (LR)]. Based on 
a 10-yr climatological assessment of 527-nm MPLNET- 
derived LR data at the Goddard Space Flight Center 
(GSFC) during May, we set LR = 50 sr to derive the ex- 
tinction coefficient. Significant variability is found within 
this climatological profile between 40 and 60 sr. This 
introduces a maximum uncertainty on solutions for the 
aerosol transmission ratio between r ref and r of 8%. 

Once [T Aaer (r re f)] 2 and \l\ A Jr)f are estimated, D 0 
estimates are independent of aerosol transmission, 
which further simplifies Eq. (11). From Eqs. (10) and 
(11), the measure normalized color ratio is 


C R (r) = 10 log 10 


( ^ 355 nora, ^ ^^355. 

I norm att ‘ 

l{^527 (0[^527 

\ norm att ; 


'ref 


'ref 


0] 2 }\ 
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where, for each wavelength, /3 Anorm t (/) is the normalized 
attenuated backscattering coefficient, defined as 


0a W 

norm att 



Px 
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PAr)r 2 
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where /l Aatt (r) = Px(r)r 2 ICf is the measured attenu- 
ated total backscatter coefficient and 6, (r) = 

2 ^ A mol att V 

/3 Aniol (r)[7\ mo] (r)] is the attenuated molecular back- 
scattering coefficient. Normalizing the attenuated total 
backscattering to the attenuated molecular backscat- 
tering (known theoretically) permits making the mea- 
surement insensitive to molecular backscattering. 

Following the Russell et al. (1979) approach and 
assuming uncorrelated variables, the measurement 
uncertainty on C R (r ) at each range gate r can be ap- 
proximated as 

[SC*(r)] 2 = [SC* (r)] 2 + [SC fi (r)] 2 , (14) 

1X 355 1X 527 


where 
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Thus, the total uncertainty in a normalized color ratio 
contains contributions from the signal measurement 
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Fig. 4. On 9 May 2012 at NASA GSFC from 0 to 5 km MSL, MPLNET (top) 527- and (bottom) 
355-nm logarithmic normalized attenuated lidar backscatter coefficient [log(Mm _1 sr -1 )] from 
1340 to 2359 UTC. 


error, the two-way aerosol transmission error, and the 
density error at both wavelengths. 

Substituting the appropriate partial derivatives, we 
have 
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Concerning the measured normalized attenuated back- 
scatter coefficient, we have 


\sp x (0 

norm att 
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where A is 355 and 527 nm. 

The uncertainty of the measured range-corrected 
signal 8[Px(r)r 2 ]l[Px(r)r 2 ] is related to Poisson photon- 
counting statistics at each range gate r (Campbell et al. 
2002). The error on lidar constant SCjf/Cjf is retrieved 
statistically with a posteriori analysis, and is about 5%. 
The uncertainty on molecular backscattering results 
essentially from uncertainties in the molecular density 
profile. In this paper, the U.S. 1976 (COESA 1976) 
standard atmosphere model is used and 3/3 Am can 
be limited to about 3% along the profile. 

The error in aerosol transmission S[r Aaer (r re f, r) 2 ]/ 
[TAaer^reG r ) : I between the range gate at the reference 
aerosol and rain-free range r ref and r, as previously stated, 
is assumed to be constant, at about 8%. The uncertainty 
on molecular transmission is again related to the mo- 
lecular density models, with S[(T Amol ) 2 ]/[(r Amol ) 2 ] « 10% 
(Russell et al. 1979). 

The effects of the aerosol contamination to back- 
scattering are discussed further below. 

3. Performance evaluation 

Light rain showers were observed late in the after- 
noon (starting at 2100 UTC) on 9 May 2012 at the 
Goddard Space Flight Center in Greenbelt, Maryland 
[38.99°N, 76.83°W; 84 m, all heights above mean sea 
level (MSL)]. Shown in Fig. 4 are the logarithms of the 
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Fig. 5. On 9 May 2012 at NASA GSFC from 0 to 5 km MSL, (top) MPLNET 355-/527-nm 
backscatter color ratio (dB) and (bottom) retrieved median droplet diameter (/xm) from 2024 


to 2248 UTC. 


normalized attenuated VIS and UV MPLNET lidar 
backscatter [Eq. (13)] composites for 1340-2359 UTC, 
respectively. Both instruments recorded data at 60-s 
temporal and 75-m vertical resolutions. The cloud base 
is shown lowering in both profiles from 5.0 to 3.0 km, 
with light rain first appearing near 2100 UTC (Fig. 4). 
The melting layer remains constant in height through 
the period; this is also confirmed by Dulles International 
Airport (38.9°N, 77.45°W; 95 m) radiosonde isotherm 
at 0°C, at about 3.1 km. A lidar dark band (e.g., Sassen 
et al. 2005; Di Girolamo et al. 2012) appears in both 
profile composites near 2.9 km, indicating a cold rain 
process as the precipitation originates from melting ice. 
Signal normalization was conducted above 3.0 km from 
1900 to 2000 UTC (60 1-min-averaged profiles). 

After 2100 UTC, precipitating fall streaks are distinct 
in both signal composites. From calculations described 
in section 2, backscatter signal profiled from precipi- 
tation at 355 nm is greater than at 527 nm. Composite 
depictions of C R and D 0 from 2024 to 2248 UTC are 
shown in Fig. 5. In the cloud layer, C R is roughly OdB, 
meaning that for cloud droplets (<30 pm in diameter) 
attenuation due to scattering is effectively the same at 
both wavelengths. The derived C R from clear sky (mo- 
lecular scattering only), being the attenuated backscat- 
tering profiles normalized to the attenuated molecular 
backscattering, is about 0 dB. 


The value of D 0 is retrieved through bivariate analy- 
sis. The range of retrieved D 0 values shown (Fig. 5, 
bottom) is consistent with light rain and drizzle, sized 
between 100 and 500 pm (e.g., Frisch et al. 1995). 

Figure 6a shows how the uncertainty relative to the 
color ratio (profile taken at 2143 UTC), calculated using 
Eqs. (14)-(18), translates into uncertainty in D 0 retrieval 
(Fig. 6b). The uncertainty in the C R is mostly dominated 
by error in lidar measurement and degrades as the 
signal-to-noise ratio. The error in D 0 is instead related to 
the unknown drop size spectrum (parameter /x) and di- 
rectly to the shape of the theoretical C R [Eq. (5); Fig. 3] 
that flattens as the raindrop diameter increases: the 
uncertainty is then bigger for larger drops (near 50% for 
600 pm. Fig. 6b). 

Errors because of multiple scattering and aerosol 
returns 

Though the 527-nm MPL instrument features a very 
narrow field of view, the UV lidar design includes one 
order of magnitude larger, and thus multiple scattering 
effects must be taken into account. O’Connor et al. 
(2004) evaluate the multiple scattering effect for a UV 
lidar with similar optical characteristics in a liquid water 
cloud consisting of relatively small droplets. Consider- 
ing some slight differences between the two UV lidars, 
multiple scattering was found to reach a maximum 15% 
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Fig. 6. (a) Measured C R profile at 2143 UTC with relative uncertainty, (b) Retrieved D 0 from C R profile and 
relative uncertainty. As the theoretical C R flattens, increasing the raindrop diameter, the uncertainty is bigger for 
larger drops, near 50% for 600 gm. 


at roughly 1 km under heavy drizzle conditions. This is an 
extreme case, and being that the attenuated backscat- 
tering coefficients in rain are lower than 30 Mm -1 sr -1 , 
a value of 10% for the multiple scattering is more likely. 
This static value has been applied for solving Eq. (12). 

At 355 nm, backscatter from aerosols is at least one 
order of magnitude smaller than that from precipitation. 
At 527 nm, the smaller backscattered power from pre- 
cipitation may be comparable to the returns from the 
aerosol particles, especially in the boundary layer. This 
will lead to D 0 being underestimated. We estimate the 
aerosol return by looking at precipitation-free areas of 
the boundary layer in Fig. 4 (top). The ratio of pre- 
cipitation to aerosol backscattered power is 5:1. For 
a measured color ratio of 4dB, the retrieved median 
equivolumetric drop size is D 0 = 140 gm. In this situa- 
tion, the aerosol backscattered power accounts for one- 
fifth of the signal and the true ratio for the rain signal 
alone would be 4.6 dB, corresponding to a true value for 
the raindrop size of D 0 = 162 gm. This value is taken 
into account to retrieve the correct D 0 in the aerosol- 
contaminated pixels mostly in the boundary layer. 

When considering each of the aforementioned sour- 
ces of algorithm uncertainty, in tandem with that due 
to signal processing, solutions for D 0 correspond with 
a maximum total uncertainty of roughly 50%. Though 


this is a considerable amount of relative error, it still 
implies that there is some reasonable measure of skill 
exhibited by the technique. Recall that, on most rainy 
days, continuous lidar profiling measurements are al- 
most immediately discarded from the negative impact 
of source pulse attenuation. Thus, the potential use of 
MPLNET measurements for estimating D Q , and with 
the thought that some uncertainty could be reduced with 
continuing instrument improvements and more stable 
calibrations, is encouraging. 

4. Conclusions 

In this paper, collocated ultraviolet and visible (355 
and 527 nm, respectively) lidar measurements are used 
to estimate the precipitation drop median equivolu- 
metric diameter D 0 over a 3-h period of light rain 
showers profiled at the NASA Micropulse Lidar Net- 
work (MPLNET)-observing site at the Goddard Space 
Flight Center in Greenbelt, Maryland, on 9 May 2012. 
Scattering calculations are performed for spherical 
liquid drops from light rain samples using Mie theory to 
link D 0 to the ratio of backscatter signals at the two 
wavelengths. These calculations are then applied to 
retrieve D 0 within precipitating fall streaks, such that 
signal attenuation effects are not limiting, including the 
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parent cloud layers to the extent of their profiling 
depth. 

We estimate that total error in solutions for D 0 is at 
about 50% for larger raindrops. The primary sources of 
uncertainty in the retrieval are attributable to the lidar 
signal measurement and the unknown shape of raindrop 
size distributions. Molecular and aerosol scattering con- 
tributions and relative transmission losses due to the 
various atmospheric constituents were evaluated to de- 
rive D 0 from the observed color ratio profiles. Multiple 
scattering is important for the UV lidar, but it exhibits 
lower overall uncertainty with respect to other error 
sources identified. 

Westbrook et al. (2010) describe the basis for the 
methodology investigated here. We find that transposing 
their technique from the near infrared to the UV-VIS 
region increases retrieval uncertainties; that is, the two 
wavelengths (355 and 527 nm) are relatively close, and 
the color ratio curve shapes become flat, especially for 
large drops (>500 /im), which increase relative error. 

Despite this, the demonstration of this method shows 
that the technique exhibits reasonable skill. With knowl- 
edge of D 0 , higher moments of the drizzle droplet distri- 
bution, including liquid water content and rainfall rate, 
can be further estimated, which will add value to the 
retrieval overall. At that stage, a value-added MPLNET 
product reporting these parameters for light rain events 
can be considered, with the goal of an improved un- 
derstanding of the microphysics and physical scattering 
properties of precipitation at near-visible wavelengths 
from global ground-based measurements. 
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